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tors, the spin and charge degrees of freedom of an electron seem to separate from 



In the strongly correlated environment of high-temperature cuprate superconduc- 

m 

each other. A similar phenomenon may be present in the strong coupling phase of 

O 

, Yang-Mills theories, where a separation between the color charge and the spin of a 



gluon could play a role in a mass gap formation. Here we study the phase structure 
of a decomposed SU(2) Yang-Mills theory in a mean field approximation, by in- 
specting quantum fluctuations in the condensate which is formed by the color charge 
component of the gluon field. Our results suggest that the decomposed theory has 

•i— I , 

an involved phase structure. In particular, there appears to be a phase which is 
quite reminiscent of the superconducting phase in cuprates. We also find evidence 
that this phase is separated from the asymptotically free theory by an intermediate 
pseudogap phase. 
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I. INTRODUCTION 

There seem to be some remarkable similarities between high-temperature cuprate super- 
conductivity in condensed matter physics and the problem of mass gap in the Yang-Mills 
theory of particle physics. It appears that in both cases the basic theoretical problem is the 
same, the absence of a natural condensate to describe the symmetry breaking that takes 
place. In high-temperature superconductors electrons do not form conventional Cooper pairs 
and the standard BCS-description of superconductivity can not be applied in any obvious 
manner. There is no obvious alternative choice of condensate that leads to superconductiv- 
ity. In a very similar way, in the case of Yang-Mills theories we do not have any natural 
candidate for a condensate of the correct dimension, that describes the mass gap of gluons. 
Could it then be that in both cases the condensate has a similar origin? 

It is definitely worth some effort to try and apply similar techniques to both problems. 
One promising method in the context of high temperature superconductivity is the slave- 
boson description, which has been studied actively This approach is based on 
the curious idea that, in the strongly correlated environment of cuprate superconductors, 
the electron (or hole) is no longer a fundamental mode of excitation, and thus electronic 
modes do not behave like a structureless fundamental object. Instead the electron can be 
interpreted as a composite particle, constructed from two quasi-particles. One of these is 
described by a charge neutral, spin-1/2 fermionic operator f ia where i is the site label and 
a =|, | is the spin index. This operator corresponds to a particle called a spinon, and 
it carries the (statistical) spin degree of freedom of the electron. The other excitation is 
described by a spinless bosonic operator bj = (frji,^)- It corresponds to a particle which 
is called a holon and it carries the electric charge of the electron. In terms of these two 
operators, the electron operator c ia decomposes as 

where we have combined the spinon operators as 

l£ • (2) 

The decomposition (JI} also introduces an internal U(l) gauge symmetry, since it is invariant 
under the simultaneous change-of-phase transformation 

h -> e%, ^ -> e i6 ^ ia . (3) 
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As a result we have a compact U(l) gauge interaction between the spinon and holon. Under 
normal circumstances we expect that the strength of this U(l) interaction increases with 
increasing energy, to the effect that at high energies the spinon and holon are confined 
into a (point-like) electron. But in a strongly correlated environment, such as in a cuprate 
superconductor, the spin and the charge of the electron can become independent excitations 
0, Q, 0] . This leads to a rather involved phase diagram, with several different regions jjj . 
One of the easiest ways to study the phase structure is using a mean-field theory. This is 
obtained by integrating over the fermions ipi a , and one finds that (c?-wave) superconductivity 
occurs when the remaining bosonic holon field foj condenses, 



Of substantial interest is also the possibility that the system can enter a pseudogap phase. 
This is a precursor to the superconducting phase with the characteristic property that even 
though the underlying symmetry is broken, the effective bosonic order parameter A& vanishes 
due to quantum fluctuations. 

Curiously, a very similar picture seems to emerge in the case of a pure four dimensional 
SU(2) Yang-Mills theory. In analogy to the slave-boson decomposition of an electron, the 
off-diagonal components of the non-abelian gluon field become composite particles, with 



Here we shall study the phase structure of the decomposed gauge theory, by following the 
mean-field approach to high-temperature superconductivity. We first construct a mean-field 
state where we integrate over the charge neutral spin degree of freedom of the off-diagonal 
gluon. We propose that in the strong coupling regime the spinless color-charge carrier of the 
gluon becomes condensed. The ensuing phase is analogous to the superconducting phase in 
cuprates. Furthermore, in analogy with cuprate superconductors we also find evidence that 
there is an intermediate pseudogap phase, a cross-over region between the superconducting- 
like phase and the asymptotically-free deconfined limit of the Yang-Mills theory. 



(b]k) = A b ± 0. 



(4) 



a separation between their color-charge and spin degrees of freedom 




II. SLAVE-BOSON DECOMPOSITION IN YANG-MILLS 



The slave-boson decomposition of the SU{2) gauge field A® (a = 1,2,3 and fi = 0, 1, 2, 3) 
proceeds as follows jj, O]: We first separate the diagonal Cartan component A 3 ^ = A^ 
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from the off-diagonal components Aj; 2 , and combine the latter into the complex field = 
+ iA 2 ^. We then introduce a complex vector field e M with 

e M e M = and e M e* = 1. 

We also introduce two spinless complex scalar fields ipi and ip2- The ensuing decomposition 
of Wf, is [j] 

W M = 4 + = M» + (5) 
This is clearly a direct analogue of Eq. ((TJ), a decomposition of into spinless bosonic 
scalars ipi,2 which describe the gluonic holons that carry the color charge of the W^, and a 
color-neutral spin-one vector which is the gluonic spinon that carries the statistical spin 
degrees of freedom of W^. 

In general, the present gluonic slave-boson decomposition is not gauge invariant. But in 
a proper gauge it can be given a gauge invariant meaning and in particular the combination 

p 2 = pi + pl = (\H 2 ) + (\H 2 ) (6) 

of the gluonic holons becomes a gauge invariant quantity. For this we introduce 

Jp 2 = f(A + pI) = / 1(a\? + (A 2 ,) 2 ] = / w,w;. (7) 

This is in general gauge dependent. But if we consider the gauge orbit extrema of (|IJ) 
with respect to the full SU(2) gauge transformations, these extrema are by construction 
gauge independent quantities. Moreover, the gauge orbit extrema of (J7J) correspond to field 
configurations W u which are subject to a background version of the maximal abelian gauge 

a, 

(d^ + igA^W^O, (8) 

which is widely used in lattice studies In the sequel we shall assume that the gauge 
fixing condition (JHJ has been implemented. The slave-boson decomposition then acquires a 
gauge invariant meaning, and in particular the condensate © is a gauge invariant quantity. 
As in (JSJ), the decomposition (jSJ) remains intact when we change phases according to 

ipi,2 -> e ie ^i,2 and e M -> e~ ie e M . (9) 

This determines an internal compact U(l) gauge structure. A compact U(l) gauge theory 
is known to be confining when the coupling is sufficiently strong [9j. The confining phase 
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is separated by a first order phase transition from the deconfined weak coupling phase. 
Furthermore, since the running of the /3-function of the compact U(l) leads to an increase 
of the coupling with increasing energy, we expect that at high energy the gluonic holon and 
spinon become confined by an increasingly strong compact U(l) interaction to the effect 
that the high energy Yang-Mills theory describes asymptotically free and pointlike gluons, 
as it should. 

But at low energy and in a strongly correlated environment, maybe in the interior of 
hadronic particles, the internal U(l) gauge interaction Q can become weak and the spin 
and the color-charge degrees of freedom of the gluon can separate from each other. If in 
analogy with (J3J the spinless color-carriers then condense 

p 2 = pI + pI = (4^i) + (4^) = a,^o , 

we have a mass gap and the theory is in a phase which is very similar to the holon conden- 
sation phase of cuprate superconductors. 

In the case of high-temperature superconductivity the basic criterion for the validity of 
the slave-boson decomposition is a dynamic one: The decomposition can occur only if the 
ensuing Hamiltonian admits a natural interpretation in terms of the decomposed variables. 
In particular, in the relevant background the holon and spinon operators should indeed 
describe proper particle states. We propose that the same criterion can also be adopted to 
Yang-Mills theories. A decomposition of the gauge field in terms of other fields leads 
to a valid description of the phase structure, only if the decomposed action has a natural 
structure and particle interpretation in terms of the new variables. In the case of Eq. (jSJ) 
this criterion turns out to be satisfied. If we write the Yang-Mills action in terms of the 
decomposed variables, it admits a natural interpretation as a two-gap abelian Higgs model 
0|. This suggests that the present Yang-Mills version of the slave-boson decomposition 
might actually identify the correct dynamical degrees of freedom that describe the non- 
perturbative phases of the theory. 

III. THE MEAN-FIELD THEORY 

In the case of cuprate superconductors, the phase structure can be investigated using 
a mean-field theory that emerges when the original theory is averaged over the electronic 
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spinon field. We now proceed in an analogous manner, and average the SU(2) Yang-Mills 
action both over the color-spinon and the Cartan component of the gauge field. Since 
we are only interested in the phase structure of the ensuing mean-field theory, it is sufficient 
to consider the free energy in a London limit where the slave-boson condensates 

pU = <hM 2 > 

are spatially uniform. 

The integration over and e M can be performed in various different ways. Our starting 
point is the one-loop result of Ref. jiol ]. which yields for the (London limit) condensates the 
dimensionally transmutated free energy 
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(10) 
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Here A is the renormalization scale and g 2 is the (A dependent) coupling constant and a 
finite renormalization A — > A sends g — > g with the familiar relation 

f = 1 , (11) 

l + f(4 ln ( A / A ) 



or in infinitesimal form 



dA ^ 3 (An) 2 ' 
The minima of (fTUj) are highly nondegenerate, and located on the p\ > 0, p<i > branch 
of the hyperbola 



\pl - plUn = A 2 exp( -^ + ^ ) • (12) 



ll£ 2 3 

Along these hyperbola, the value of the free energy is 



11 .4 ( 48vr 2 22 \ 

-k exp -tt-j + — • (13) 



2416tt 2 "V H(7 2 3 

Since (II Oj) is a one-loop approximation it can not be used for providing numerically 
accurate predictions. For this, high-precision Monte Carlo simulations are needed. But if 
one is only interested in the qualitative features of the phase diagram, the explicit form (fTUj) 
is adequate, as shown below. 

Indeed, if we assume that the Yang-Mills /5-function has no zeroes so that 

dx 

W) <0 • 
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the minimum values (|12|) and (|13|) can be represented in the renormalization group invariant 
form 



9 

2,2 _ \ i ' ' dx 22 



(Pi - P2)Ln = A 4 exp( -4 / — + - 



and 

9 9 9 

F min = -4A 4 • y • exp ^-4 J + = -A J (3(x)dx ■ (pi - p 2 2 ) 2 n 



mm' 



Consequently we expect that the qualitative features of our conclusions have a validity which 
extends beyond the one-loop level. For the present purposes it is sufficient to start from the 
notationally simpler version 

1 



f = ^(pI-pD 2 -[i + x 



In — a 

A 2 



(14) 



We normalize F with the factor (1 — Act), set A 2 = 1 and redefine 

A 1 * 



21- \ a 

and arrive at the final version of the free energy that we shall use in our analysis: 

F = \{pl-pl?-(i + \-Hp\-pl)^ . (is) 

In figure ^ we have plotted this free energy for A = 100, on the entire (pi,p2) plane. The 
generic features of this potential, a ridge along the lines p\ = ±p 2 , and a narrow hyperbolic 
valley on both sides of these lines, are independent of A, but the depth of the valleys and 
steepness of the potential are more prominent for larger values of A, as used here. 
The Landau pole at 

\ = l/a (16) 

separates the strong coupling region with A > Ija from the weak coupling region with 
A < 1/a. Since the latter region includes the small coupling limit of the original Yang-Mills 
theory and since the free energy (fTUj) can only be reliable for weak coupling, we shall in 
the following concentrate on the region < A < 1/a. Notice that in terms of the redefined 
coupling A in ()15|) . this corresponds to the region of a positive A. In particular, if a in ()14j) is 
large, the strong coupling limit A — > oo in (|15|). does not necessarily correspond to a strong 
coupling limit of the original model (|14|). 
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FIG. 1: The free energy F in equation (|15[). for A = 100. The physical branch corresponds to the 
quadrant pi^ > 0. [color online.] 

IV. CLASSICAL ASPECTS 

We first consider the properties of the free energy (JT5j) at a classical level, where we do 
not include the quantum fluctuations in the spatially uniform London limit condensates pi 2- 
This free energy has the following classical scaling symmetry 

P 

P = ~ > 

c 

A' A 



l + 2Aln(c 2 ) ' 

F = ~ A ' l + 2Alnc 2 ' ^ 
We can employ this scaling symmetry to restore the parameter A in the free energy; see 
(|14)1. Indeed, it is obvious that this scaling symmetry reflects the renormalization 

group symmetry of the original Yang-Mills theory, with the scaling transformation of the 
coupling constant A a version of (JTTJ). 

In addition, as a function on the entire (pi, p 2 ) plane the free energy has a discrete 
symmetry since it only depends depends on a polynomial combination of the condensates 

F(px,P2) = f(( P l-pl) 2 ). 
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and in addition we also have the gauge invariant polynomial combination in Eq. (J7J), 

p 2 =pl + pl . (18) 

We are interested in linear transformations that act on the (pi, p 2 ) plane and leave both 
polynomials intact. These polynomials are exactly the basic invariants that generate the 
octonic dihedral group D4 (also called 4mm), which is the nonabelian symmetry group of 
the square. 

The four branches of the hyperbola that minimize (|15j) . 

(Pi - Pi? = exp ( , (19) 



A 

are separated by (non-analytic) ridges along the lines p\ = ±p 2 , and mapped to each other 
by the transformations. At the minima along the hyperbolic valleys the free energy is 
given by 

E min = -lexpK^J . (20) 

This ground state is highly degenerate, but the combination on the left-hand side of (|T9*|) 
is not the proper gauge invariant condensate. The gauge invariant condensate is given by 
Eq. ()18j) . and we can employ it to remove the infinite degeneracy of the hyperbolic vacuum: 
From (|19J) we conclude that the ground state value p 2 = v 2 of the gauge invariant con- 
densate (|18|) is bounded from below by a non-vanishing quantity, 

P 2 = Pi + p\ = v 2 > \p\ -p\\= exp f i— -J . (21) 



2A 

When v 2 is larger than the lower bound in (J21j) . there are eight solutions (pi, p 2 ) to the 
equations that define the vacuum 

2,2 2 

Pi + P2 = V > 

pi-p\ = ±exp (L=±) . (22) 



2A 

But when v 2 coincides with the lower bound there are only four solutions, 

pi — ±v & p 2 = , 
p x = & p 2 = ±v . (23) 

which correspond to the vertices of the hyperbola. The solutions are mapped onto each 
other by the dihedral group D4, and selecting any one as the ground state breaks the -D4 
symmetry. 
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The solutions of (|22|) describe the generic situation where both condensates are non- 
vanishing. The solutions are /^-degenerate, but we remove this degeneracy when we select 
the (physical) px,2 > quadrant. The remaining ground state is doubly degenerate under 
exchange of p\ and P2, which correspond to the physical scenario that in general the London 
limit densities are unequal. 

Finally, the degenerate solutions ()23|) correspond to the limit where one of the two con- 
densates vanishes, and again by selecting the physical quadrant p 12 > we remove the 
degeneracy. 

According to (|2"T| the ground state value of (fT%|) is non-vanishing for all non-vanishing 
values of the coupling constant A. This suggests that in the Yang- Mills theory the gauge 
invariant condensate (JBJ) is also nonvanishing for all non-vanishing values of the coupling 
constant. This would mean that the mass gap in the Yang-Mills theory is present for all 
values of the coupling, and it vanishes only asymptotically in the short distance limit where 
the gluons become asymptotically free and massless. 

V. QUANTUM MECHANICS - NUMERICAL APPROACH 

The classical treatment of the mean-field theory in the previous section suggests that 
the condensate (jBJ) is always non-vanishing, hence a mass gap is present for all nontrivial 
values of the coupling. We now want to inspect what effects spatially homogeneous quantum 
fluctuations around the classical mean-field value have on this condensate. For this we need 
to improve the free energy so that it also includes the contribution from the momenta 
7ii 2 that are canonically conjugate to the (spatially homogeneous) condensates p\^- For 
computational simplicity we consider these condensates to be defined over the entire (pi, p^) 
plane. This results in a -D4 symmetry, and by selecting the physically relevant values p\^ > 
for the condensates we then break this discrete symmetry. 

The conjugate momenta are the generators of spatially homogeneous translations. Their 
inertia is undefined, and we therefore add a parameter M. The improved free energy can 
be interpreted as a Hamiltonian 

h = ^(*1 + A) + \(p\-pI) 2 (i + ^(pI-pI) 2 ) • (24) 
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It corresponds to the effective action 

T 



S c s = J dt(n 1 p\ + ix 2 p2 - H[ir, p\) , (25) 
o 

and the equations of motion for S c s are invariant under the following extension of the 
scaling transformation (|17|) 

/ P 
P = ~ , 
c 



1 7T 



7T 



A' 



1 + 2A ln(c 2 ) c 2 
A 



l + 2Aln(c 2 ) ' 
i! = v / TT2Ahi(c 2 ) • ct . (26) 

We also note that the action is clearly invariant under the dihedral symmetry group. 

In order to study the effects of quantum fluctuations in the condensates we investigate 
the solutions of Schrodinger equation 

~ (91 + dl) HPuP2) + \{p\- PIT (l + A In (p? - pif) ^ Pl ,p 2 ) = E^(p u p 2 ). (27) 

We have studied this Schrodinger equation (j2T|) numerically, using a highly-accurate finite 
difference approximation, on grids of varying size and spacing, using up to 400 x 400 grid 
points. We have analyzed both the ground-state wave function and several of the low- 
lying excited-state wave functions when the coupling constant A in (|27p varies for fixed M. 
According to the relation between ()14j) and (|15j). this surveys the phase structure of the 
theory at couplings below the Landau pole. 

We note that since the Schrodinger equation (J27|) is invariant under the action of the 
dihedral -D4, the wave functions can be chosen to have definite -D4 transformation properties. 
Unfortunately, most discussions of point groups, see e.g. fl2| . look for representations in 3D 
space, where one can distinguish between the groups C± v and D4, but these groups act 
identically in the xy plane. Following Mulliken's (A, B ID irrep; E 2D irrep) or Koster's 
notation (Yj) as discussed in Ref. [l^, we have five possible representations of this group in 
two dimensions, see table E 

In that table we have also listed representative wave functions for all the irreps. Looking 
at the symmetry of the wave functions, we expect Ax, Bi and one of the E cases with zeroes 
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TABLE I: Character table (trace over classes of elements of the representation matrices) for D±. 
The five classes on the top line are, respectively: E identity; 2C^ Z rotations over ±7r/2 around the 
z-axis; C 2z a rotation over tt around the z-axis; 2a y reflections in the x or y-axis; 2a xy reflections 
in the lines x = ±y. 

E 2C^ Z Ciz 2a y 2a xy representative wave function 
Ax = T 1 1 111 1 f(x 2 + y 2 ) 
A 2 = T 2 1 1 1-1-1 xy(x 2 - y 2 ) 
B l = T 3 1 -1 1 1-1 x 2 -y 2 
B 2 = T 4 1 -1 1-1 1 xy 

E = T 5 2 -2 (x, y) and (x + y, x - y) 

on the lines x = ±y to form four almost degenerate states as A grows large, as borne out by 
figure El below. 

When A — > the Schrddinger equation ()27|) reduces to the x 2 y 2 model which has been 
studied in detail in [1 3l | . In particular, it has been established that the spectrum of the x 2 y 2 
model is discrete, the eigenstates are normalizable, and the ground state energy is separated 
from E = by a non- vanishing gap. 

In figure El we depict the behavior of the numerically constructed ground state wave 
function for different values of A for M = 1. Very similar behavior is found for other values 
of M, but as analyzed in more detail below, the similarity is greatest if we compare solutions 
for identical values of AM. We find that the wave function exhibits three different kinds of 
qualitative behavior. There is a weak coupling region < AM < A a ~ 10, an intermediate 
coupling region A a < AM < X b w 500 and a strong coupling region X b < A. In all cases the 
ground-state wave function lies in the lowest symmetric representation (A^ of D 4 . These 
regions have the following characteristic features: 

a. Weak coupling: In the weak coupling region with MA < A a ~ 10, we find that the 
ground state wave function is qualitatively reminiscent of the ground state wave function 
in the x 2 y 2 model, in particular it has a single maximum which is located at the origin of 
the (pi, P2) plane. We also find that the value of the wave function at its maximum varies 
very slowly as a function of A, especially for a large value of M, which means a more tightly 
localized wave function, see figure (JH1); Both the shape of the ground state wave function and 
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X=50 X=\00 X=200 



FIG. 2: The ground state wave function for A from to 200. A result for a larger value of A is 
given in figure |SJ [color online.] 

the location of its maximum suggest, that in this weak coupling region quantum fluctuations 
tend to restore the system towards the symmetric state p 2 ~ so that there would not be 
any mass gap in the underlying Yang-Mills theory. 

Clearly, we find non-zero condensate values for any finite value of M. As suggested by 
figure EJ the condensate goes to zero at A = 0, as M goes to infinity but there remains a 
cross-over to a broken phase at stronger coupling. This is consistent with the fact, that in 
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FIG. 3: The upper panel shows the two condensates p 2 and {{p\ — p 2 .) 2 ) 1 / 2 , and the distance of 
the maximum value of the ground state wave function to the origin. The lower panel gives the 
value of the ground state wave function at the origin. The solid lines show the values for M = 1, 
and the dashed lines for the larger value M = 8. [color online.] 

the limit of vanishing coupling the Yang-Mills theory describes free massless gluons. 

Independent of M, when A approaches A a the value of p\ overshoots the classical value 
given by the right-hand side of (flUjl . Consequently the expectation value (P) detects the 
presence of symmetry breaking and the ensuing nontriviality of the condensate, even though 
this is not reflected in the location of the maximum value of the ground state wave function. 

Such a behavior where the condensate © detects a symmetry breaking while the wave 
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FIG. 4: The condensate ©, as a function of MX for various values of M, as indicated in the plot, 
[color online.] 

function tends to retain the symmetry, is reminiscent of the pseudogap phase P . As a con- 
sequence we propose that in the weak coupling region < A < A a the underlying Yang-Mills 
theory is in a pseudogap phase, a cross-over region which terminates in the asymptotically 
free theory as A — > 0. Presumably this region is intimately related to a Coulomb-like phase 
in the Yang-Mills theory. 

b. Intermediate coupling: When 10 ~ MX a < MX < MXb ~ 500 there is a clear qual- 
itative change in the behavior of both the ground state wave function and the condensates 
(JHJ). For these values of the coupling the origin in the (pi, Pz) plane becomes a local mini- 
mum, and instead there are now four maxima in the wave function. These maxima are all 
located the same distance p max from the origin, and related to each other by the -D4 symme- 
try. Both the value of the ground state wave function at the origin, and the condensate (jOJ 
decrease essentially linearly in the logarithmic scales of figure El while the value of p max very 
rapidly approaches the classical limiting value determined by (fHJj) . as A — > A&. In particular, 
as A — > Xb the value of the condensate © becomes less than its classical bound in (fT§)). but 
is still clearly bounded from below. 

In this intermediate coupling region both the ground state wave function and the conden- 
sate (jBJ) behave similarly, and in a manner which suggests that the underlying Yang-Mills 
theory has a mass gap. Indeed, the behavior is quite reminiscent of the superconducting 
phase in cuprate superconductors. We find it natural to propose that this region of the cou- 
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pling constant describes a superconducting mass-gap phase of the Yang-Mills theory, maybe 
a magnetic dual to the confinement phase. 

c. Strong coupling: When A — > A& ~ 1000 we detect a new transition, towards a strong 
coupling regime At < A. Now the value of the ground state wave function essentially vanishes 
at the origin, see figureEl The value of the condensate (jOJ again increases, and asymptotically 
approaches the value p max = exp(— 1/4) which is the classical A — > oo lower bound value 
(12 lj) for the minimum distance between the potential minimum and the origin. Indeed, for 
the entire strong coupling region A& < A we find that the difference between the classical 
and quantum values of the condensate is very small, suggesting that in this region one of 
the condensates essentially vanishes. Consequently as A — > oo the system becomes driven 
towards a degenerate ground state where one of the condensates asymptotically vanishes, 
while the other becomes asymptotically determined by the classical theory. 

The strong coupling region retains the major characteristics of the intermediate coupling 
region: There is a mass gap, and the ground state wave function is peaked at a nontrivial 
value of the condensate, even though in the infinite coupling limit one of the (quantum) 
condensates seems to vanishes asymptotically- it seems that the two quantities p 2 and 
(0°i ~ pV) 2 ) 1 ^ 2 coincide as A — > oo. But in this region the ground state wave function has the 
additional characteristic property that it (essentially) vanishes in a neighborhood around 
the origin, thus becoming (essentially) separated into four disjoint components. This means 
that the lowest four states, consisting of two one-dimensional representations and one two- 
dimensional one, become degenerate. While we do recognize that in a finite dimensional 
quantum mechanical model there always remains a (vanishingly small) tail of the wave 
function at the origin, the numerically observed vanishing of the wave function in the vicinity 
of the origin is very definite. Consequently we envision, that in the underlying field theory 
with its infinite number of degrees of freedom, there is a true transition where the tunneling 
between the four dihedrally symmetric branches of the ground state wave function becomes 
totally suppressed. 

Finally, we find that the at these high values of A the wave functions of the four lowest 
states become degenerate, see figures El The symmetries of these states correspond to A X) 
E and B\ irreps, as argued above. 
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FIG. 5: The four lowest eigenstates (from left to right and top to bottom) of the problem (|27|) for 
A = 400. [color online.] 

VI. QUANTUM MECHANICS - ASYMPTOTIC ANALYSIS 

The present model is a generalization of the x 2 y 2 model, a notoriously complex system 
|l^ |. While the A —>■ limit of our numerical results reproduce the known properties of the 
x 2 y 2 model, there is a need to confirm the main features of our A / results by formal 
analysis. For this, we now consider the relevant asymptotic behavior of the ground state 
wave function. 



A. Large distance behavior of the wave function 

For A = the behavior of the solutions of (071) are known and have been discussed in 
: The spectrum is discrete, the eigenstates are normalizable, and the ground 
state energy is separated from E = by a non-vanishing gap. Our numerical investigations 
suggest that these conclusions persist for non-vanishing values of A. We now proceed to 
verify this using asymptotic analysis of the Schrodinger equation (j77|) . In particular, we 
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wish to establish that the wave function is indeed normalizable. 

When pi ^ ±p 2 the potential in (|27j) is bounded from below by a positive quadratic form. 
Consequently any peculiar, unexpected behavior in the ground state wave function must be 
concentrated near the lines where p\ = p\. This is best studied in hyperbolic coordinates 

Q 

£=~{p\-pl), and v = PiP2 ■ (28) 

We note that even though these coordinates only cover half of the (pi, p 2 ) plane, we can use 
them to study the full behavior of the wave function for large values of pi ~ ±p 2 - In these 
coordinates the Schrodinger operator is 

H = -Ve + V*(dx + d m ) + 2^ 2 (1 + Aln(4e 2 )) 

and we are particularly interested in the behavior of the wave function for large values of rj 
and small values of £. 

We first consider the known A = case, this leads us to the Schrodinger equation 

We wish to implement an asymptotic separation of variables. For this we select rj to be large 
and positive (alternatively large and negative). Since the potential depends on £ alone, we 
can introduce a transformation in this variable that allows separation to (almost) take place: 
We write 



and 

With this Ansatz we get 

yg(y) (-/" (*) + 2x 2 f(x)) - 



x = £/ri 1/4 & y = V 1/2 



<t>{£,v) = f{x)g{y) 



f(x)g"(y) | l 

4 y 



5xg(y)f'(x) f(x)g'{y) xf(x)g'(y) 



16 4 4 

x 2 g(y)f"(x) 



16y 2 

Thus, the problem is separable to leading order in 1/y, 



Ef(x)g(y) . (29) 



f"(x)+2x 2 f(x) = rf(x) (30) 
+ ryg(y) = (31) 
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The function / is clearly one of the harmonic oscillator states, and for the ground state of our 
Schrodinger equation we must have the lowest energy eigenstate of the harmonic oscillator. 
The equation for g then becomes 

g"(y) 



4 



+ V2yg(y) = 



since for large values of y the value of the energy becomes irrelevant. For a normalizable 
wave function, the only acceptable solution is 

g(y) = M{2^y) (32) 

which decays rapidly. This is consistent with our numerical simulations, and the (A = 0) 
results in |l3j: The "tendrils" of the wave function along the potential valleys are indeed 
decaying very rapidly. In original coordinates, 

0( Pl , p 2 ) » exp(-i(p? - Pl) 2 /(pip 2 )) Ai(2 5 / 6 ( PlP2 ) 2 ) (33) 

We now consider the general case: As above, our approach is based on asymptotic sepa- 
ration of variables obtained by rescaling £, which allows us to combine the terms multiplying 
/" with the rescaled potential. Let us therefore look at the £ —>■ at; scaling of the general 
potential, as studied in Eq. (JTTj), 

y A (0 = 2£ 2 (l + Aln4£ 2 ) . (34) 

We find 

V x (aO = a 2 (l + X\na 2 )2e (l + + A In 4^ = , (35) 

where 

P = a 2 (l + Alna 2 ) & A = A - . (36) 

(1 + a X) 

If we now make a dependent on rj, and define x = £/a(?y), we get a matching condition by 
requiring a common ^-dependent factor for the leading second derivative with respect to x 
and the potential 

rja 2 = a 2 (l + A In a 2 ) 

which has the solution 

v 1/4 

2r] \ 
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FIG. 6: The functions a (left) and (3 (right) vs. A and [i. [color online.] 

where plog is the "product logarithm" (inverse to xe x ). 

The functions a(rj) and (3(r]) (see Fig. are increasing functions of rj, and will provide 
us with a second expansion parameter. Now we separate variable in x and r\ as before, 

<P(£,v) = f( x )g(v) , 

which leads us to the following generalization of (j2ll|) 
- V9"(v)f(x) + (3(rj)9(v) (-f"(x) + V~ x{v) (x)f(x)) - 

m^nf g{T]) {3xf{x)+X 2 fix)) 

- P(vW{v) x f'( x )g'(v) + 7;P(v)a"(v)xf\ x )g(v) = Ef(x)g(r]) . (38) 

We now ignore all but the first three terms- it can be verified numerically that all other 
terms are small-and separate variables 

- f"(x) + V~J{x) = e(X(v))f(x) , (39) 

-vg"(v) + P(v)<Kv))g(v) = o . (40) 

The value /3(r/) is larger than r/ 1//2 , its value when A = 0, as can be seen in Fig. El The 
eigenvalue e(A) is larger than \/2, if we take rj large enough so that A = 0(1), see Fig. [7| 
Since A — * for rj — > oo, we then have a rapid decay of the wave function for large 7], as 
expected. 

The equation for f(x) is interesting for other reasons as well; as shown in figure |S1 we 
find that for large values of A the wave function separates into two parts. Remembering that 
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FIG. 7: The lowest eigenvalue of the ID Schrodinger equation (|39|) as a function of A. 
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FIG. 8: The lowest eigenfunction of (|39|) as a function of A. [color online.] 

x = corresponds to the lines pi = ±p 2 , this supports our assertion that the wave function 
separates in 4 disjoint parts for large A. 



B. Large values A > A^ 

In the region of large coupling, the hyperbolic valleys of the free energy become very deep. 
We are interested in the asymptotics of the ground state wave function, when it becomes 
separated into four disjoint components. We continue to utilize the hyperbolic coordinates, 
but we shall now expand around the minima of the free energy. 

In hyperbolic coordinates, the minimum of the free energy 
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remains very close to £ = (less than ^ = 0.303 . . . which is the value for A — > oo). Con- 
sequently the previous asymptotic analysis remains valid, and we can immediately conclude 
that the wave functions are decaying rapidly. 

However, it is also of interest to consider the limit of a very deep potential directly, and 
for this we introduce coordinates from the minimum, scaled with £ , 

£ = £o(:r±l) & r 1 = ^y . 

We then expand in powers of x, keeping leading terms only. The Hamiltonian simplifies to 

H = ~^V(T+xY+?(d xx + d yy ) + 2\e (x + l) 2 [log(x + l) 2 - 1] 

^-^TTf(d xx + d yy )+2\e (-l + 2x 2 ) . (41) 



£ M 

With (j, = A£q and e = E^ + fj, we this leads to the eigenvalue problem 

- ^ (d xx <j) + d yy <P) + 4/ix 2 = 60 . (42) 

We now wish to consider the properties of solutions to this Schrodinger equation: We sub- 
stitute 

<P{x,y) = f{z)g{y) , 

with 

z = x(l + y 2 ) 1/8 ■ 

When we ignore terms containing lower order or mixed derivatives in addition of powers of 
x, we find that the equation takes the form 

(1 + y 2 Y /4 9(^d zz f - Wf) + f^TT7^d yy g = efg . 

Assuming again the lowest harmonic oscillator eigenstate for /, 

/ = exp(-(z/6 2 ) 2 /2) , 



with 



we find for g, 



b z = (2/iM)" 1 / 4 , 

^^d yy9 + V2jl(i + y 2 ) 1/ '9 = eg . 
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In order to obtain an analytic solution we shall assume fi to be so large that we again can 
make the harmonic approximation. This gives 

-d yy g + \fu T j2y 2 g = (e - y/2^)g . 

Thus 

g = exp(-(y/b y ) 2 /2) , 
with by = (8/M/i) 1 / 8 , and the ground state energy is given by 

e = (2/i/M) 1 / 2 + (/i/(2M) 3 ) 1/4 . 

The energy for the original problem can thus be expressed as 

(2 / i/M)V 2 + ( / i/(2M) 3 )V 4 - / . 

1 ' 

So 

which is a good approximation only when the wave function fg has no overlap with those 
from the remaining three valleys - and when our harmonic approximations are valid. 

We also conclude that the scaling of the condensates observed in the previous section is 
indeed taking place; the width of both / and g depends on the combination /iM = AM£q, 
and for large A £ is approximately constant, leading to the observed scaling in AM. Since 
the wave function contract slowly to the maximum point (most slowly for g), we find that 
indeed we have in the limit A — > oo the case where one of the two condensates disappears, 
as stated above. A further numerical analysis using the separated wave function confirms 
that this approach is very slow, and we must go to extremely high values of A to see the 
point where we can't distinguish between the maximum and the expectation value of p 2 . 

VII. CONCLUSIONS 

In conclusion, we have investigated the phase structure of pure SU(2) Yang-Mills the- 
ory using a slave-boson decomposition of the gauge field. We have employed a mean-field 
approximation where we account only for spatially homogeneous fluctuations in the gluonic 
holon fields. Our analysis suggests that the decomposed theory has an involved phase dia- 
gram, resembling that of cuprate superconductors. At intermediate couplings, there seems 
to be a gapped phase which is separated from the asymptotically free high energy limit 
by a pseudogap phase. Furthermore, we find that a mass gap appears to persists in the 
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strong coupling limit even though asymptotically one of the two holon condensates appears 
to vanishes. 

In our analysis we have employed a version of the maximal abelian gauge. In this gauge we 
have the advantage, that many results are available from first principle lattice simulations; 
see for a review. In particular, it has been observed that the (electric) confinement of 
color relates to the condensation of magnetic monopoles in the dual Higgs phase. Here we 
have inspected a (renormalization group invariant) perturbative one-loop approximation to 
the Yang-Mills effective action, in terms of decomposed variables that have a natural mag- 
netic interpretation. Our results do not account for topologically nontrivial configurations, 
consequently it is not directly clear how our results could relate to the (electric) color con- 
finement as observed in the lattice simulations. For a comparison, we need a first principles 
numerical lattice analysis in terms of the separated spin and charge variables. We also need 
a better understanding of electric-magnetic duality in terms of these variables. However, 
even at the level of the (crude) approximation that we employ here, it appears that when 
formulated in terms of the separate spin and charge variables the Yang-Mills theory has a 
very rich phase structure, not easily described in terms of the conventional gluonic variables. 
Our results suggests that the possibility of a spin-charge separation in the Yang-Mills the- 
ory may occur, and deserves to be addressed by extensive first-principles lattice simulations. 
Furthermore, there is a need to address theoretical issues such as electric-magnetic duality 
and the description of the Yang-Mills theory in terms of the spin-charge separated dual 
variables. 

Indeed, if gluons can become decomposed into their independent holon and spinon com- 
ponents, it could have deep consequences to our understanding of the fundamental structure 
of matter. 
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